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This article provides a convenient framework for quantitative evaluation of the entanglement gen- 
erated when two structureless, distinguishable particles scatter non-relativistically in one dimension. 
It explores how three factors determine the amount of entanglement generated: the momentum dis- 
tributions of the incoming particles, their masses, and the interaction potential. Two important 
scales emerge, one set by the kinematics and one set by the dynamics. This method also provides 
two approximate analytic formulas useful for numerical evaluation of entanglement and reveals an 
interesting connection between purity, linear coordinate transformations, and momentum uncertain- 
ties. 
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Two structureless, distinguishable particles scatter non-relativistically in one dimension. Specifically, the two parti- 
cles are prepared independently, interact via a potential with finite range, and separate again. When particles with no 
initial classical or quantum correlations scatter, they generally emerge from the interaction region entangled. In this 
context we ask the question: how much entanglement between the particles is generated and on what does it depend? 
The topic of entanglement in particle scattering has intrinsic interest: scattering is a fundamental dynamical process 
and entanglement is the fundamental resource for quantum information processes. Despite the potentially complicated 
dependence of entanglement generation on the momentum distributions of the incoming particles, their masses, and 
the interaction potential, this article identifies some general features of entanglement generated in scattering systems. 

Entanglement generation in non-relativistic scattering of structureless, distinguishable particles has been considered 
previously [H, [3, 0, 0, H, H, 0, B Hi ■ Most previous treatments consider interactions of a single species [2, 4, 5, 6, 7,"^. 
Some consider scattering by particles confined in a harmonic well interacting by delta- function potentials [2,, _6, 9] and 
were directly inspired by experiments with ultracold atoms in traps. Analytic results for the energy levels of such 
systems [13, Ell can be used as a discrete basis for calculating entanglement. Others have considered Gaussian wave 
functions of unbounded particles interacting via specific limited-range potential models, including hard core J_, 
"cavity" (hard core plus delta- function) and double delta- function Q. Each of these have considered different 
aspects and have particular strengths. This research seeks to unify some of the different observations of these papers 
and provide a useful framework for analyzing general properties of entanglement generation in scattering. 

Scattering also supplies a rich context for exploring techniques of quantum information theory with continuous 
variable systems. A bipartite continuous variable pure state will have a wave function f{xi,X2)- The state is 
unentangled if the wave function is separable, and then it can be written as f{xi, X2) = g{xi)h{x2)- An entanglement 
measure will quantify the non-separability. Some inseparable wave functions can be written as 



where each of the /; have disjoint domains of support and are therefore orthogonal. Each fi can be thought of as the 
wave function of a "mode" of the system, and entanglement due to the superposition of modes can be called "mode- 
like" [l^]. Either the modes fi are separable, and then the only entanglement is mode- like, or they are themselves 
irreducibly inseparable. Use the purity of the one-particle reduced density matrix pi2 defined below, and the purity 
of wave function like H]) is 
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In other words, the purity of a state that can be reduced into orthogonal modes is the sum of the purities of those 
modes. 
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Using these ideas, several different analytic expressions for the purity will be given for different approximate or 
incomplete-knowledge situations that arise in particle scattering. Scattering can split the wave function into modes, 
such as a transmission mode and a reflection mode, or it can distort the wave function shape. Either effect can render 
it inseparable. In resonance scattering > the rapid variation of the scattering amplitude with relative momentum 
leads to such distortion, but even when the amplitudes are constant, the wave function can be rendered inseparable 
by the distortion caused by reflection. Entanglement that originates from kinematic effects, like entanglement due to 
reflection, can be distinguished from dynamic effects using the methods described here. 

As an example, consider hard core repulsion in one-dimension. In 4], the case of two particles with the same mass 
is studied. The amount of entanglement is measured by the inverse of the purity of the one-particle reduced density 
matrix. Its variation with time is numerically calculated, and first it grows from no entanglement to a maximum 
value as the particles approach and interact. Then it recedes back to no entanglement. In contrast, as first shown in 
[l[ and also discussed in [1, Q, if the particles have different masses, then the outgoing states will remain entangled 
after scattering, and possibly very entangled, unless 



where ai — Api is the momentum uncertainty [l6[. In one-dimensional scattering hard core repulsion means only 
reflection can occur, and this leads to the following interpretation: During the time interval the scattering actually 
occurs, even with equal masses, the two-particle wave function is distorted by interference between the parts of the 
wave function that are still incoming and the reflected parts. For the equal mass case, after a long time the final 
wave function will return to an undistorted separable state. However, if the masses are different, the process of 
reflection will engender correlations between the particles that remain after the particles leave the scattering region. 
Reflection means reversing the internal momentum, but different parts of the single-particle momenta distributions 
will be affected in different ways by this transformation and so the wave function will be distorted. More details are 
given below. 

While the approach of this article uses the S-matrix, thus excluding exploration of the time dependence of the 
entanglement, it does allow the separation of kinematic and dynamic entanglement effects by exploiting symmetries 
and conservation laws. Another critical feature of our approach to entanglement in one-dimensional scattering is the 
application of the mathematics of tensor product structures (TPSs). The notion of entanglement depends on how 
a system is factored into subsystems. Most treatments of entanglement in scattering are concerned exclusively with 
interparticle entanglement, in other words entanglement with respect to the TPS 



where Hi is the Hilbert space of a single particle. Other TPSs are useful for scattering entanglement. For example, 
one can define a TPS which factors the into the internal and external degrees of freedom (see [12, 13]). In [l^] it was 
proved that for non-relativistic scattering, entanglement with respect to the internal-external TPS is unaffected by 
scattering because Galilean invariance interactions will factor into unitary operators operating on each product space. 
Generally, a change of observables leads to a new notion of entanglement, and this can occasionally be exploited to 
achieve useful analytic and approximate results. In particular, the Schulman condition ^ can be derived as a special 
case of the purity of Gaussian wave functions under linear transformations of the observables, and the entangling 
distortion of reflection can be analytically calculated in the same framework. 

The outline of this paper is as follows: First, we describe the how to quantify entanglement for two particles in one 
dimension using the purity of the reduced density matrix, and we show how coordinate transformations lead to new 
TPSs and new definitions of entanglement. Then we describe the scattering dynamics and show how the boundary 
conditions of scattering imply that the total entanglement of the scattered state is partially mode-like, in the sense 
that it is the sum of the the entanglement of the transmitted and reflected states. Next we consider several different 
analytic approximations to the purity and the conditions under which they apply. In the process, three different 
aspects of dynamical entanglement will be thrown into distinction, entanglement due to superposition of transmission 
and reflection, entanglement due to wave packet distortion under reflection, and entanglement due to wave packet 
distortion due to convolution with a rapidly-varying scattering amplitude. At the end, we will comment on possible 
applications and indicate directions for further research. 
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II. QUANTIFYING ENTANGLEMENT OF TWO ONE-DIMENSIONAL PARTICLES 



For a system of two spinless, structureless, non-relativistic particles in one dimension, a complete set of commuting 
observables (CSCO) is the momentum operators for each particle {Pi, A}- The generalized eigenvectors (eigenkets) 
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of these operators are |pi,P2) and they have delta- function "normaUzation" 

(Pi,P2bi,P2) = ^(Pi ~Pi)S{p2 -v'-i)- (5) 

A general state of the two-particle system 1$) G 7i can be realized by its two-particle momentum wave function 
0(Pi,P2): 

\(j)} = J dpidp2(t)ipi,P2)\pi,P2}- (6) 

The CSCO used above implies a tensor product structure, i.e. a way of partitioning the total Hilbert space 

n = Hi(E)H2--L^{R)(E}L^{R)=L^iM'^), (7) 

where i^(R) is the space of Lebesgue square-integrable functions In turn, this tensor product structure (TPS) 
implies a notion of entanglement: interparticle entanglement. 

To quantify the interparticle entanglement, we use the purity of the one-particle reduced density matrix pi2, which 
for the sake of brevity we refer to as the interparticle purity: 

Pi2(</')-Tri(p2) = Tr2(p2) (8) 

The one-particle reduced density matrix is 

pi = Tr2(p) = J dpidp[dp2C^{pi,P2W{p'i,P2)\pi){p'i\. (9) 
Inserting ^ into ([5]) gives a simple, symmetric expression for the interparticle purity: 

Pl2((/>) = / dpidp2dp[dp'2(l){j)l,P2)<P*{p'l,P2)<pij>'l,p'2W{Pl,p'2)- (10) 



Assuming that the pure state \(j)) is normalized, the purity will be in the set pi2 G (0,1]. For a separable state 
4's{pi,P2) = 0i(pi)02(p2), one calculates pi2 = 1: 



Pl2{(t>s) = I dpidp2M'^P2'/'l(Pl)</'2(P2)'/'t(pi)02(P2)0l(yi)02(y2)0l(Pl)</'2(P2) 

dpi\Mpi)A ( f dp2\MP2)\A ■ (11) 



One could imagine calculating the entropy of entanglement E12 = — Trpi log pi for a pure continuous- variable state, but 
diagonalizing the one-particle reduced density matrix is less amenable to producing analytic results. To diagonalize, 
either one calculates the eigenvalues of the continuous reduced density matrix, which is typically a hard problem, or 
one discretizes the momentum spectrum in some way, which could obfuscate some of the very effects we are trying to 
uncover. 

Other CSCOs are compatible with interparticle entanglement. For example, the CSCO of the particle position 
operators {Xi, X2}, leads to the same TPS ([7]) as the CSCO {A, ^2}- The CSCOs and (therefore the wave functions) 
are related by Fourier transformations. These are local unitary operators with respect to the interparticle TPS 
(O. The interparticle entanglement of a position wave function (l){xi,X2) can calculated by Fourier transforming the 
position wave function to the momentum variables and the evaluating the purity pi2, or equivalently by evaluating 

Pi2{(t>) ^ / dxidx2dx[dx'24>{xi,X2)<p*ix[,X2)(t>{x[,x'2)(t>*{xi,x'2). (12) 



Similarly, global translations in position space and global boosts in momentum space are local unitary operators 
with respect to the interparticle TPS. These coordinate transformations will be represented on ([7]) by non-entangling 
operators because the generators of these transformations factor into a direct sum of local operators. For example, 
a global translation in space by a is generated by the total momentum operator P = Fi <^ Ii + I2 €5 A ■ When 
this is exponentiated, the translation operator factors into local unitaries: exp(iaPi) ^ exp{iap2) [isj . Similarly, one 
could boost the system without generating interparticle entanglement, including a global boost to the frame where 
the expectation value for the total momentum P is zero, i.e. the center-of-mass (COM) reference frame. However, 
the interaction Hamiltonian does not separate into a sum of local operators on the interparticle TPS ([7]) since the 
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interaction potential depends on the relative positions. As a result, time translation, and therefore the scattering 
operator, are not local operators and can change the amount of interparticle entanglement. 

There are other CSCOs which are useful for understanding the scattering of two particles in one dimension, but 
which are not compatible with the TPS ([7]). For example, the COM and relative position coordinates 

Xc — ^Ji-iXi + ^ji2X2 (13a) 

Xr — X2 — Xi 

and conjugate momentum coordinates 

P = P1+P2 (13b) 
where fii — mi/(mi +^2), are the natural coordinates for most scattering calculations. The COM-relative momentum 

csco 

{p = A ® ii + 12 ® A, Q = /i2A ®h- Mill ® A} (14) 

defines the internal-external (IE) TPS for one-dimensional particles 

n = Hc(^nr L'^iM)^ L^{R) = L'^iR^), (15) 

and the COM-relative position observables {Xc, Xr} define an equivalent TPS. 

Note that the total Hilbert space in ([7]) and ([15]) is the same; only the partition into subsystems is different. The 
subsystem partitions ([7]) and (|15p are isomorphic, but they lead to inequivalent definitions of entanglement. The 
interparticle wave function can be transformed into the IE wave function by 



(t>iP,q) = / dpidp2(t)ipi,P2){p,q\pi,P2) 



dpidp2(j){pi,P2)5{pi - ^J.lP - q)5{p2 - f^2P + q) 

= Hf^iP + q,f^2p- q)- (16) 

As can be seen from the above expression, a pure state \<j>) that is separable with respect to the single-particle TPS 
([7]) may not be separable with respect to the IE TPS and vice versa. 

The IE TPS is useful because the Hamiltonian for all Galilean invariant interactions will have the form 

H = P^/{2M) + Q^/{2m) + V{Xr) 

= Ac ® Ii + Ii ® Hr (17) 

where M = mi+m2 and m = 17111712/ M = 111H2M. This separation of variables is exactly why this coordinate system 
is used to solve the Schrodinger equation, but it also can be seen as a statement about the dynamical entanglement 
of scattering. Since the Hamiltonian is a sum of local operators with respect to the IE TPS fTS]) . the entanglement 
between the relative and COM degrees of freedom (or between the internal and the external degrees of freedom) is a 
dynamical invariant for any Galilean invariant interaction. This is a result that also holds for any Galilean invariant 
interaction in three dimensions and for particles that have internal structure such as spin, although the definition of 
the internal and external CSCO changes [13, [l^l • 

The transformation of the momentum variables (jl3bp is just one example of a real linear transformation of the CSCO. 
Using the notation p = {pi,P2) to denote a vector in (pi, p2)-space, we consider how to calculate the entanglement 
of a given wave function (j){pi,p2) = 0(p) with respect to new variables z = Tp, where T is a real, 2x2 matrix 
with detT = ±1 to preserve normalization. These variables are the eigenvalues of a CSCO {Zi,Z2} given by the 
transformed momentum operators Z = TP. This CSCO in turn implies a new TPS 

and therefore to a new notion of entanglement is measured by a transformed purity. The purity with respect to the 
new TPS is 

PT12(<A) = dzidZ2dz[dZ2(f){zi,Z2)(f>*{z[,Z2)(t>{z[,Z2)(f)*{zi,z'2), (19) 
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where 0(2:1,22) = "^z)i,(T ^7)2) and z = (zi, 22)- 

As an example, consider a state that is the product of two Gaussian distributions of momentum: 

(j)a{p) ^ NiN2e''''''^e ^"'i e^^^^^e "'2 , (20) 

where Ni = {2T:af)~^^^ , ki are the central momenta, are the central positions, and ui are the momentum uncer- 
tainties for each particle's Gaussian. Because it was constructed as a product of single particle wave functions, the 
interparticle purity of this state pi2 — 1. However, the transformed purity (|19p of \(I)g) for a transformation matrix 

T=f;:V (21) 



where \ru — st\ — 1, is 



Pti2((/>g) = '^^'^^ — (22) 

See the Appendix for the derivation. The quantity \22\ does not disappear even when momentum the variances Ci 
are small. The scale for this kind of entanglement is set by the ratio of the momentum uncertainties, as can be seen 
more clearly by defining c — 1T2 /cri : 

Pti2('/>g) = J. „ „ „ „ „ :■ (23) 
^(r^ + s^c^){t^ + v?c^) 



Specifying to the matrix that transforms from the particle momentum coordinates to the IE CSCO, 



1 1 

/U2 -/Wi 



(24) 



we find the IE entanglement is 



Pie{<Pg) = „ „ „ =■ (25) 

Inspecting this function, one finds that a separable Gaussian state will have entanglement with respect to the IE TPS 
unless the condition 

^ = ^ (26) 

is satisfied. This is equivalent to the Schulman condition ([3]) which was found by explicitly considering the transformed 
wave function [l| 



Pl ^2 



(p-fei-fea)' 



(PG 

xe V 1 2/ e V''"! '^"■2] (27) 

and looking for the separability conditions. The relevance of this fact will be seen when the kinematics and dynamics 
of scattering are discussed below. 

III. DYNAMICAL ENTANGLEMENT THROUGH SCATTERING 

The in-state \4>in) of a scattering system is the state of the particles a long time before the interaction, formally in 
the limit t ^ —00. The boundary conditions of scattering put certain constraints on the functional form of 0in(pi,P2), 
and these affect the entanglement that can be generated. 

First, it is assumed that the two particles are prepared independently and that before the scattering they have no 
correlations (classical or quantum). Therefore, the pure in-state \(j)in) is separable with respect to the interparticle 
TPS (HI and the wave function can be factored 

|0m) = |0l) ® |'/'2) ^ 0m(Pl,P2) = 0l(pi)0(p2) (28) 
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as a product of pure, single-particle states G Tii. Since a global boost in momentum space is a local unitary 
operator, we can perform our analysis in the reference frame where the expectation value of the COM momentum P 
is zero. Therefore, without any loss of generality, we will assume that the expectation values of the single particle 
momentum operators in the state 0i„ are (Pi) = k and (P2) — ~k. Further, we will also assume that the single- 
particle wave function (/)i{pi) has support only (or to a good approximation) in the range pi > and that 02 (P2) has 
support when p2 < 0. Then the particles are moving toward each other, and eventually transmission and reflection 
can be treated as modes. This condition is implies 

dpMp)<p2ip) = 0. (29) 

As a final assumption, the particles are initially located far from each other, or if the potential has a finite range, 
outside that range. An example of such a state that satisfies these criteria is the product of Gaussian momentum 
wave functions (|20p with central momenta fci = — fc2 ~ k a positive number, initial positions a2 = —ai = a large and 
positive and where ai/k ^ 1: 

(f>G{p) = NiN2e'^P^-P'^''e e , (30) 

When (Ji/k = 1/6, the integral in ((29|) is less than a millionth. 

In elastic scattering, the out-state {(pout), formally the state in the limit t 00, will also be an element of Hi ®'H2- 
The scattering operator 5* is a unitary operator on 7i that connects the in-state and out-state 

\(t'out) = S\(t)in). (31) 

If the scattering potential Vr depends only of the relative positions and has finite range (and also for more general 
assumptions) , one can show that the scattering operator exists and its properties can be derived from the Schrodinger 
equation [l^. While S preserves the purity of the total state, it is not a local operator with respect to the TPS The 
interaction will generate correlations between the momenta and we general expect that (j)out{pi,P2) 7^ 4''i{Pi)'t>'2{P'2) ■ 
Then the essential problem becomes, for a given S and what is the interparticle purity pi2{4>out)- 

In one dimension, the S-matrix (the realization of the scattering operator in a certain basis) has a particularly 
simple form. Using the IE CSCO, it is a sum of a transmission term and a reflection term: 

{p,q\S\p' ,q') = S{p,p',q,q') = Se{p,p')S^{q,q') 

^ S{p' ~p){tiq)6{q~q')+r{q)S{q + q')), (32) 

where we have factored it into an external (p-dependent) part and an internal (q-dependent) part to emphasize that 
it is a local operator in the IE TPS. The functions t(q) and r(q) are the transmission and reflection amplitudes and 
by unitarity |i(9)P + |9(9)P = 1- These can be determined by solving the 1-D Schrodinger equation in the relative 
momentum coordinate using the reduced mass; analytical forms exist for a variety of potentials. Using this expression 
for the S-matrix, the out-state can be written as 



I'/'out) = / dpdqdp'dq'S{p,p,q,q)(t>inip,q)\p,q) 

dpdq {t{q)(j>,r,{p, q) + r{q)(f>i,i{p, -q)) \p, q) 

dpdq (0tra(p, q) + (prefiP, ?)) \p, q) 

dpdqcpout {p, q) |p, q) (33) 
To calculate the interparticle purity of the out-state, one must convert (j)outip,q) into 4iout{pi,P2)- 

|0o«t) = J dpidp2(j)outiPl,P2)\Pl,P2} 

= J dpidp2\pi,P2){<PtraiPl,P2)+(l}ref{Pl,P2))\Pl,P2)- (34) 

and then insert this into PH)) . One finds that 

4>traiPl,P2) = t{qi2)(t)in{Pl,P2) (35) 
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where qi2 = ^2Pi - A*iP2 and 

(l)ref{Pl,P2) = r{ql2)(p^n{Pl,P2) (36) 

where Pi = (/^i — ^i2)pi + 2/xip2, and P2 — 2^2Pi + (M2 ^ f^i)P2- The expressions for arise due to the effect 
of transforming the particle coordinates to the COM-relative momentum coordinates, flipping the direction of the 
relative momentum, and then returning to the particle coordinates. In terms of matrices, this transformation looks 
like 

p = Mp = T-iFT„nP (37) 
where F is a matrix that flips the relative momentum. In components, this is 



M 



fJ-2 -I J \ -I J \fl2 

{pi - M2) 2/^1 

2^2 {P2 - fJ-l) 



(38) 



One can prove that the transformation psp maps disjoint regions of the (pi,p2)-plane into disjoint regions. In 
particular, it reflects the line pi — —p2 about the origin without change of scale, i.e. {k,—k) — > {—k,k). This 
means that 4itra{pi,P2) and 0re/(pi,P2) have disjoint domains of support; i.e. they are 'modes' in the sense of 
the introduction. When calculating the integral in pi2 (0out), one takes the four- fold product of (pout with different 
permutations of pi, p[, p2 and pg- The only terms that survive the integrals over all four variables are those terms 
which involve only (ptra or only (pref- In other words, the integral for pi2{4>out) breaks into two parts, the purity of 
the transmitted state and the purity of the reflected state: 

Pui'Pout) = PuiPtra) + Pui'Pref)- (39) 

Note that since the probability for transmission T = WcptraW^ and reflection R — ||'/>re/|P are less then one, one can 
show pi2{(j>tra) G (0, T^] and Pi2[(pref) G (0, i?^], so pi2{(j>out) is bounded from above by one. 



IV. USEFUL APPROXIMATIONS 



Very few assumptions went into (j39p : the interaction was Galilean invariant and sufficiently well-behaved that 
the S-operator exists, and the in-state wave function must satisfy certain boundary conditions typical of scattering 
experiments. If pin and S are explicitly known, then pop can be calculated numerically. However, there are further 
useful analytic approximations that are applicable under reasonable circumstances. 



A. Entanglement of reflection 

The first approximation applies when the transition amplitudes t{q) and r{q) do not vary much with q. Note 
that since we calculate in the COM frame, the transformation (jl3bp implies the expectation value of g is fc and the 
uncertainty of the relative momentum q is Ag = Y^(/xiApi)2 -f (/i2Ap2)^- If we assume that t{k) ^ Aqdt/dq, then 
we can approximate the functions t{qi2) and r{qi2) by their values at the expectation value of relative momentum, 
tik) and r[k). With this assumption, the purity of the transmitted state can be calculated immediately: 



Pl2{Ptra) = j dpidp2dp[dp'2ptraiPl,P2)Ptra{Pl:P2)Ptraip'l,p'2)(PtraiPl:P2) 

= |i(fc)|^ J dpidp2dp[dp'2pin{pi,P2)P*n{p'i,P2)Pin{p'i,p'2)PiniPl^P'2) 

= \t{k)fpl2{P^n)^T^ (40) 

because puiPin) = 1. If t{q) does vary rapidly, then the integral cannot be simplified because the factors like t{qi2) 
will generally lead to a marked distortion of the wave function that does not preserve separability. This is particularly 
true when there is a resonance and t{q) varies to a maximum value of 1 over a small range of g 



8 



The purity of the reflected state cannot be simphfied in such a way even if r{q) is a constant. The inversion of the 
internal momentum q leads to entanglement between pi and P2- The reflected wave function is 

<PrefiPi,P2) = r-(gi2)(/)™((Mp)i, (Mp)2) = r-(gi2)(/)^(pi,p2), (41) 

where M is the matrix ([55)1 that enacts the transformation on {pi,p2) due to q ^ —q. One can think of this as another 
example of a linear coordinate transformation described in Section II, except now this is an active transformation of 
the wave function instead of a passive transformation of the coordinates. One can use the formula (fT9|) to calculate 
the covariance properties interparticle purity for a state 4^{pi,P2)- By making the the substitution dpidp2 — dpidp2, 
one finds 



Pl2('fe) =mi2(0««)- (42) 

Assuming a Gaussian in-state of the form ([30|) , one then finds 

This is a function of the masses and momentum uncertainties of the in-state single-particle wave functions, but it can 
also be re-expressed entirely in terms of the ratios 1122/ mi and (T2/<^i- The scale of the entanglement by reflection is 
set by the relative masses and momentum uncertainties of the scattering particles. The purity takes its maximum 
value (and therefore there is the least entanglement) when the masses are either equal mi = TO2 or when the Schulman 
condition mi/af = m2/cr| is satisfied. In [8], the authors identify these same two conditions as necessary for non- 
entangling scattering from hard core potentials. Since hard core scattering assumes only refiection for all q, we see 
that any entanglement generated by hard core scattering is of a purely kinematic origin and independent of the energy 
scale of the dynamics. Reflection is an entangling process except for equal mass scattering or for when the Schulman 
condition is satisfied. This effect does not disappear even in the limit of narrow momentum distributions. 

For all other potentials, there can be both transmission and reflection, and therefore both kinematic and dynamic 
entanglement even in the constant amplitude approximation. The total interparticle entanglement in this approxima- 
tion can be efficiently calculated using the sum 

Pl2{^out)=T^ + R^Pl2i^). (44) 

with T + i? = 1 and pi2(<%^) given by (I43p . Note that the strength of the dynamical entanglement depends on the 
transition probabilities, which in turn depend on the energy scale of set by the scattering particles and the interaction, 
unlike the kinematic entanglement of reflection which is independent of the overall energy scale. 



B. Coarsest approximation: Two two- level systems 



In any experiment, there is incomplete measurement of the out-state and one does not have access to full information 
about the wave function. Real experiments cannot measure probability distributions over continuous variables; the 
continuous variable must be broken into 'bins' with some resolution. This coarse-graining is equivalent to discretizing 
the system. In the simplest model where the continuous variable of interest is broken into TV disjoint bins, then each 
one-particle system may be approximated as an iV-dimensional discrete system. Then the minimum purity possible 
for a coarse-grained bipartite continuous- variable system would be whereas the 'unbinned' purity is limited from 
below by zero. 

For example, in the worst case, it might only be possible to detect whether after scattering the distinguishable 
particles remain on their original sides or they have transmitted past each other. Then one cannot measure the all 
correlations in the two-particle wave function 4'out{PiiP2) that may exist, but one only knows that the function has 
support around one of two domains: the transmitted mode centered around (fc, —k) and the reflected mode centered 
around (—k,k). Then the particles are effectively two level systems. In this approximation, each particle has two 
states that we label as |±), the sign of its momentum in the COM frame. For the in-state, we have 

= 1+)®!-) = I + -). (45) 

The out-state will in general be a superposition of the transmitted and reflected states 



|0o«t) i\ + -) +r| - +). 



(46) 
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The numbers t and f are the effective amplitudes chosen such that \t\'^ = T and |fp = R. If the transmission and 
reflection amplitudes are relatively constant over q around fc, then t « t{k) and f « r{k), but this approximation 
works in the more general case as long as the (f>ref{pi,P2) and 4'tra{pi,P2) are essentially orthogonal. 

In this approximation, the space of states for each particle is effectively realized by and the total Hilbert space 
n for the system is reahzed as 7^ = Hi ® x = C^. We can calculate the interparticle entanglement by 

looking at the reduced density matrix for one of the particles: 

Pi =Tr2(|<^)(0|). (47) 

For the in-state, we have 

P1=(J2) (48) 

and then we can calculate the entanglement by the purity of the reduced density matrix pi2 = Tti{pi) = 1. As 
expected there is no entanglement. For the out-state, we have 

Pi^(j^). (49) 

The entanglement of the particles in this "no wave function" approximation is 

Pi2=T^ + R^- (50) 

Interparticle entanglement is maximized when T ~ R = 1/2. Then pi2 — 1/2, which is the same maximum en- 
tanglement as for two two-level systems. Physically, the superposition of transmission and reflection leads to the 
entanglement of the particles, and when these components have equal likelihood, the particles are most entangled by 
the collision. Not just correlation is required for interparticle entanglement (the in-state had perfect anti-correlation 
between the particle velocities), but indeterminacy between the two modes of transmission and reflection is required 
if one cannot measure the correlations between the particles within a mode. 

V. CONCLUSION 

Our primary concern has been the interparticle entanglement of the pure out-state. However, there are several 
factors that limit how measurable this entanglement could be. First, typically the in-state will be a mixed state pm- 
The boundary conditions on the in-state still require it to be separable, but purity will no longer be an acceptable 
measure of the entanglement generated in the out-state SpoutS^ ■ Second, scattering experiments extract information 
about the out-state from transition rates and cross sections. The measuring devices can be represented as operators 
A, and then measurements are connected to the probabilities by the Born rule prob — Tr{pout^)- Since there are no 
perfect measuring devices, the kind of coarse-graining described above will always reduce the measurable entanglement 
from its theoretical maximum. Finally, as extensive research on Bell inequalities has shown, some of the correlations 
measured in a scattering experiment are classical, or they may be quantum but do not rise above some threshold that 
could not be explained by a local, realistic theory. See Q for further discussion of this last point. 

Aside from issues of direct measurement of entanglement due to scattering, another interesting avenue of research 
is the connection between entanglement in individual scattering events and entanglement properties of systems of 
particles. For example, scattering dynamics give a contribution to the equation of state for effective many-body 
Hamiltonians of systems like ultracold atomic gases. Such analyses usually assume that the exact form of the interac- 
tion is unimportant, but as seen above, kinematic effects may still drive the generation of entanglement in scattering. 
For example, L.S. Schulman discusses in |3|] the case of two gases of atoms with different masses interacting through 
scattering. As a result of the kinematic entanglement due to reflection, he derives the result that the interparticle 
entanglement for any interspecies pair reduces to zero after multiple scatterings and that the particle wave functions 
tend toward Gaussian distributions that satisfy the Schulman condition ([3]). How to implement this hypothesis ex- 
perimentally is unclear, but advances in the control of ultracold atomic systems promise future novel experiments in 
interparticle entanglement. 
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Appendix 

Here we present a derivation of 122]) . Our method is similar to that found in Q which is inspired by squeezed-state 
quantum optics calculations, but we use Gaussian wave functions, not Gaussian Wigner functions. A useful result 
will be that for an N x N diagonal matrix B with elements and an TV-dimensional real vector x, we have the 
Gaussian integral 

/ d~xexp(-x^Bx) = nJl^- (51) 

J — CO V " 



First we rewrite (l20t as 



(p) =7ViA^2exp(za-p)exp(-i(p-kfV(p-k)) (52) 



where Ni — {2'Kaf) a = (ai, a2), k = (fci, and V is the 2x2 matrix 





2o-^ / 

The new coordinates are defined by Tp = z. These coordinates imply a different CSCO {Zi, Z2} and therewith a 
new TPS. We want to calculate the entanglement of ([5^ in this new TPS. Using the notation T = T^^, the wave 
function in the new CSCO is 

0(z) = iViiV2Cxp(ia-fz)exp(-i(tz-k)^V(fz-k)) 

= iViiV2exp(ia-fz)exp(-iz^Vz)cxp(-kVtz)cxp(-ik^Vk), (54) 

where V = f^VT. 

As a next step, we define the following variables: z = (21,^2), z' = (zJ,Z2), z = (zJ,Z2), z' = (zi^z!^). Then, the 
purity Pti2 is 



Pti2(0g) = J rf'zd^z>(z)0(z) cj){z')<P{z') 

= N^N^ j d2zd2z'exp(ia-t(z-z + z'-z')exp(-2k^Vk) 

X exp(-i(z'^Vz + z^Vz) + (z')^Vz' + (z')'^Vz') exp(-kVt(z + z + z' + z') 



N^N^ j Szd?i' exp(-2k^Vk - 2kVf (z + z') 



X exp(-i(z^Vz + z^Vz + (z')^Vz' + (z')^Vz'). (55) 
To integrate this in the fashion of ([5T|l . we define the following variables: w = T — k, y = z — z', and 



V' 



22 



where Vn = u^Vn +t^V22 and V22 — s^Vn +r'^V22- By defining D as just the diagonal elements of V, the exponential 
terms in ()55p can be re-expressed as two purely quadratic terms: 



PtM^g) = \NtNl j d2wd2yexp(-i(w^Vw))exp(-l(y^Dy)), (57) 
where the fact that 4(i^zd^z' — J (PwcPy has been used. Finally, using ([5T|) . pxi2 can be evaluated 



, , . 1 ,r4.r4 /27r 2tt 2tt 2tt 



(58) 



V'(rV2-fs2(T2)(iV2+u2a2^ 
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Note that in [J position uncertainties are used instead of the momentum uncertainties and so their a is the inverse 

of our a. 

[17] Actually, as is normal in quantum mechanics, we will only consider elements of a dense subset $ C 7i whose elements 
are realized by wave functions that are well-behaved, i.e. smooth, rapidly-decreasing. An example is the Schwartz space. 
This restrictions allows the use of the Dirac basis vector expansion and the use of Riemann integrals instead of Lebesgue 
integrals in calculations. 

[18] A local translation like exp(iaiPi) dg) exp{ia2P2) is also non-entangling, but this operator does not commute with the 
interaction Hamiltonian and it is not a generator of the Calilean symmetry group for the whole system. 



